Chiral Spin Liquid in Correlated Topological Insulator 
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In this paper, we investigate the topological Hubbard model - the spinful Haldane model with on- 
site interaction on honeycomb lattice with spin rotation symmetry by using slave-rotor approach and 
find that chiral spin liquid exists in such a correlated electron system of the intermediate coupling 
region. By considering the anyon nature of excitations, chiral spin liquid may be the ground state of 
the topological Hubbard model. The low energy physics is basically determined by its Chern-Simons 
gauge theory. 
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The Fermi liquid based view of the electronic proper- 
ties has been very successful as a basis for understanding 
the physics of conventional solids including metals and 
(band) insulators. For the band insulators, due to the 
energy gap, the charge degree of freedoms are frozen. If 
there exist spontaneous spin rotation symmetry break- 
ing, the elementary excitations are the gapless spin wave 
and the gapped quasi-particle (an electron or a hole) that 
carry both spin and charge quantum numbers. However, 
in some special insulators, the elementary excitations 
with fractional quantum numbers of an electron may ex- 
ist. People call them quantum spin liquid states[l-3]. 
There exist different types of ansatz of spin liquid: Z 2 , 
U(l), SU(2) and SU(2) x SU(2)[2, 3]. These different 
spin liquid states have the exactly the same global sym- 
metry, as conflicts to Landau's theory, in which two states 
with the same symmetry belong to the same phase. In 
particular, there exist quantum spin liquid states break- 
ing time reversal symmetry, of which the elementary ex- 
citations are anyons with fractional statistics. People call 
them chiral spin liquid (CSL)[3, 4]. There are two types 
of CSLs - the abelian CSL with abelian anyonic excita- 
tions and non-Abclian CSL with non- Abelian anyons. 

Recently, nonAbelian CSL state has been predicted 
in the Kitaev model on honeycomb lattice or in its 
generalizations [5]. On the contrary, although Abelian 
CSL has been proposed much earlier than non- Abelian 
CSL, till now people don't know any types of model with 
the (abelian) chiral spin liquid as the ground state. Then 
one issue here is may people realize chiral spin liquid in 
certain many-body systems? To answer above question 
we study the quantum properties of the so-called topolog- 
ical Hubbard model on honeycomb lattice with spin rota- 
tion symmetry by using slave-rotor approach and propose 
that chiral spin liquid may exist in such a correlated elec- 
tron system of the intermediate coupling region, of which 
there exist anyonic excitations. 

The topological Hubbard model on honeycomb lattice : 
The Hamiltonian of the topological Hubbard model on 
honeycomb lattice is given by 



H = Hv 
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Here -Hr- is the spinful Haldane model as 
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t and t' are the real hopping between the first-neighbor 
and the second-neighbor on the different and the same 
sublattices, respectively. e ZVij is a complex phase into the 
second-neighbor hopping, and we set the direction of the 
positive phase is clockwise | = \ ) [6]. U is the on-site 
Coulomb repulsion, a are the spin-indices representing 
spin-up (a —\) and spin-down (a =1) for electrons, fi 
is the chemical potential and fx = [7/2 when the system 
is half-filling (in this paper we only study the case of 
half- filling) . and denote two sites on a first- 

neighbor and a second-neighbor link, respectively, 
and foil are the number operators of electrons with up- 
spin and down-spin respectively. It is obvious that there 
is spin rotation symmetry but no time reversal symmetry 
for the topological Hubbard model. 

For free fermions (the on-site Coulomb repulsion U 



is zero), the spectrum is Ek = ±y (tjy) + where 



|£k| = ty3 + 2cos(v / 3fc y ) + 4cos(3fe x /2)cos(v / 3fc !/ /2) 
and = It' ^2 sm (k • bj) . The parameters ai , a.2 and a3 

i 

are the displacement from one site to its nearest neigh- 
borhoods and bi = a2 — a3, h?, = a.3 — ai, etc. The 
length of the hexagon side has been chosen to be unit. 
One can see that there exist an energy gap at 
the points ki = ^L(l, and k 2 = ^(-1, -775). Due 



to the existence of nonzero TKNN number [7], there ex- 
ists the integer quantum Hall effect a xy = =j-. Here the 
parameter 2 comes from the contributions of electrons of 
up-spin and down-spin. Therefore, for the free fermions, 
the ground state is a topological insulator with quantized 
anomalous Hall effect [6, 8, 9]. 

An issue is whether the topological insulator is stable 
for the interaction case. To examine stability of the topo- 
logical insulator against on-site interaction, we will use 
the slave-rotor approach to study the topological Hub- 
bard model. Slave-rotor approach has been widely ap- 
plied to study the quantum liquid states near Mott tran- 
sition of correlated electron systems [10- 15]. By the slave- 
rotor approach, we find that chiral spin liquid appears of 
the intermediate coupling region. 
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FIG. 1: (color online) Phase diagram at T = . There are 
four regions : I is TI state, II is the quantum spin liquid, III is 
an AF order with QAH effect, IV is the trivial AF order. The 
blue line and black line are (¥) c i and (t") C 2, respectively. 



Slave-rotor approach : In slave-rotor approach, we rep- 
resent electronic operator Cj a into Cj a — e l6j fj a where 
the fermion spinon fj a represents the spin degree of free- 
dom and e l9j represents the charge degree of freedom 
together with slave- rotor's constraint ^2 fjufja = 1- 

Here we introduce an additional variable - the angular 
momentum L — ihdg associated with a quantum 0(2) 
rotor 6. Then the Hamiltonian in Eq.(l) turns into 

Heff = "* E (flfjrXlXj + h.C.) - flJi* ( 2 ) 
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where hi is a Lagrange multiplier for slave-rotor's con- 
straint, pi is a complex Lagrange multiplier for \Xi\ =1 
(X i = e^). 

We introduce four variational parameters Q f = 

Eflfja) , Q' f = 

a I ij{nn) 

.Ee^flfj.) ■ To 
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obtain the five parameters Qx, Q/, Q'x, Q'fi P, we solve 
the following equations self-consistently, 
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FIG. 2: (color online) (a) : The energy gaps of rotor Ax 
(the red line) and fermionic spinon A/ (the blue line) of the 
topological Hubbard model at T = and t' = 0.12. (U/t) cl 
and (U/t) c2 are the quantum phase transitions, (b) The non- 
zero spin chiral order parameter. 



where 



3k 



4 cos (3k x /2) cos(V3k y /2) + 2 cos(v / 3fc y ) (4) 



and e k = -Q x |£k| 
of unit cells. 



t'Q' x gk. N s denoting the number 



After the calculation, we find that the topological in- 
sulator is stable below a critical interaction strength, 
f < (t)c1 (See the results in Fig.l) With increasing 
interaction strength, we get non-zero solutions of Qx, 
Qf, Q' x , Qj, p. As a result, the ground state turns into a 
quantum spin liquid state characterized by a finite gap of 
rotor excitation. The excitations are not fermions gener- 
ated by c|, instead, they are rotors and fermionic spinons. 



Then we get the mean field effective Hamiltonian as 



H, 



eff 



Hf + Hx where 



h, = -t e QfflM - 1' E Q'f^flha + h. 



c. 



(5) 
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and 



H x = -tJ2 QxX*X j -t> ]T Q'xXfXj + 

(6) 

+ / l ^i l + p^|X 4 | 2 + / l .c. 



After diagonalizing the Hamiltonian in Eq.(5), we get 
the energy spectrum for the two-flavor fermionic spinons 

E f = ±^Q}\^\ 2 + Qf\^\ 2 . One may get the energy 

gap for the fermionic spinons as Af = 6V3t'Q' f . On 
the other hand, the spectrum of the charge excitations 
Ex is Ex = + U Ek. In the quantum spin liquid 

state, Ax is not zero, A.y = (p + min (ffk)); when 

approaching the phase transition between topological in- 
sulator and quantum spin liquid state it becomes zero 
due to rotor condensation. From Fig. 2. (a), we can notice 
that the rotor's gap is much larger than the fermionic 
spinon's gap (the energy scale of Ax is 10 times to that 
of A f ). 

By this method we get a quantum spin liquid state 
with fermionic spinons by adding the on-site interaction 
to the spinful Haldane model. When further increas- 
ing the interaction strength, the quantum spin liquid 
is unstable against antifcrromagnctic (AF) spin density 
wave (SDW) order. Such AF-SDW order is described by 
(4,oA°-) = |(1 + (-l)VM). Here M is the staggered 
magnetization. In HF mean field approach, we may get 
the self-consistency equation for M by minimizing the 
ground state energy. To characterize different orders of 
the topological Hubbard model (the topological insula- 
tor, the quantum spin liquid, the AF-SDW), we plot a 
phase diagram in Fig.l. ( y) c2 denotes another critical in- 
teraction strength that divides the quantum spin liquid 
and the AF-SDW (See Fig.l). In particular, in Fig.l, 
one may see that there exists a narrow window between 
quantum spin liquid state and the trivial AF order - an 
AF-SDW order with quantized anomalous Hall effect. 

Effective Chern-Simons theory : In the following parts, 
we will focus on the quantum spin liquids between topo- 
logical insulator and AF-SDW state (region II in Fig.l). 

In the quantum spin liquid state, because the rotor 
excitation Xi has a big energy gap, we may integrate out 
it and concentrate only on the spinon excitations. The 
fluctuations of hi and the phase fluctuations of Qf, 
amount to coupling the fermionic spinons to a compact 
U(l) gauge field by the minimal prescription. After 
considering the fluctuations around the mean field saddle 
point, we get the effective model of fermionic spinons 



with U(l) gauge invariance 

L f = fja(dr - ia T ,j +h - n)f ja (7) 
-tQ f ]T e^flf ja 

(id)," 

-t'Q' f ]T e^e^fjj^ + h.c. 

where /io = (hi) = /i. Hence the continuum version of 
above model becomes the two flavor massive Schwinger 
model with the Lagrangian as Cf = itf)^ ^(d ^ — ia^ip + 
rmpip where m = Af/2 is a fermion mass and ip a = 
i^alo = ( ha A, ffaB, UaA, fiaB ) and a = 1, 2 labels 
the two points k x = ^(1, -^) and k 2 = ^(-1, -^). 

Considering the quantum fluctuations of fermionic 
spinons, we get a two dimensional dynamics Maxwell 
model of the gauge field a M as C a = -^(d^a^) 2 where 

e 2 , = The compact U(l) gauge theory is always 

confining [20]. However, the induced Chern-Simons (CS) 
term will lead to deconfmement. Integrate over fermions 
by using 1/m (gradient) expansion approach, we obtain 
the CS term C cs = YW\^ e ^ Xa ^ ax where N = 4 
(two- flavor plus two spin components) [16, 17]. 

Finally we obtain an effective CS theory of the quan- 
tum spin liquid state with the Lagrangian C c g = Cf + 
C a + Ccs- 

Chiral spin liquid : After obtaining the effective CS 
theory, the quantum spin liquid state ( region II in Fig.l) 
is identified to chiral spin liquid. Such a topologically 
ordered spin liquid breaks time-reversal symmetry while 
preserves all other symmetries (spin rotation symmetry, 
translation symmetry, ...). 

Firstly we point out that the quasi-particle is really 
anyon. Let us consider single 7r-flux excitation ($ = tt) in 
the quantum spin liquid as shown in Fig. 3. The vacuum 
expectation value of the fermion number (N { ) is related 
to the spectral asymmetry of the Dirac Hamiltonian 

1 f 00 1 1 

(N l ) =-- alE -Im Tr (- ——) sign(E) (8) 

2 J_ OQ tt Hf-E-ie 

where Hf is the Hamiltonian of the fermion spinon [18]. 
And the fermionic number is also related to the Atiah- 
Patodi-Singer invariant r\ H = — i(A f ) which represents 
the difference between the number of states with posi- 
tive and negative energy The Atiah-Patodi-Singcr index 
theorem states that due to the quantum anomaly the 
fermionic number of the Dirac operator, equals the topo- 
logical charge as (N l ) = — ^-S-^r, N = 4. It represents 
a fact that a 7r-flux excitation with half topological charge 
carries one fermion number, |(iV')| = 1. That means a 
7r-flux excitation is really a bound state of 7r-flux and 
a fermionic spinon. Due to nontrivial AB phases upon 
adiabatic exchange of charge and flux, 7r-flux turns into 
scmion - special type of Abelian anyon. By the fusion 
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FIG. 3: (color online) Scheme of an anyon as a composite 
object of a fermionic spinon and a 7r-flux. 1, 2, 3 correspond 
three vertices of a equilateral triangle to define the spin chiral 
order parameter in Eq.(9). 

rules of Abelian anyons, one may find a statistical angle 
is | = f . 

Secondly, we calculate the topological degeneracy of 
CSL, a topologically ordered spin liquid. In the tempo- 
ral gauge, ao = 0, and on a torus, the ground states are 
characterized by zero momentum gauge fields (a x ,k=o, 
a Vi k = o). After straightforwardly calculations [3], we may 
get the effective Hamiltonian of (a x ,k=Oi Oyk=o) a s 

nj = CPe*-Ag y ) 2 (Vey-AeJ 2 _ A _ a B ,k=o 

ne ff — 2A/ X + 2M y wnere j\ 6x - 2n , 
A dy = ^ and M x = M y = (L x and L y 

are the lengths of the system along x- and y-directions, 
respectively). This model corresponds to a particle on 
a plane with a finite "magnetic field". The strength of 
the "effective magnetic field" is obtained as B e ff = Jj- 
There exists a two-unit flux tube through the center of 
the torus. So the degeneracy is given as V = 2. 

Next we calculate the edge states of the CSL from the 
effective CS theory. We know that the charges of a M are 
quantized as integers. Then the effective CS theory has 
two right-moving edge excitations. The two branches of 
the edge excitations are described by the following ID 
fermion theory£ odgc = Y^a^aR^t ~ v R d x )ip aRi where 
a = 1,2. iJ} a R carries a unit of a M charge. That means we 
get spin-charge separated edge states : the edge modes 



carry only spin current [21, 22]. 

Thirdly, an important property of CSL is the non-zero 
spin chiral order parameter which has a non zero expec- 
tation value in a phase with broken P and T symmetry. 
Spin chiral order parameter is a rotationally invariant 
operator defined through [4, 19] 

X(123> = (Si • (S 2 x S 3 )> (9) 

If the sites 1, 2, 3 correspond three vertices of a equilat- 
eral triangle in a plaquette, we may estimate the mean 
field value of x ( i 23> to be i (sin$) \Q' X \ 3 where <I> is 
the gauge invariant flux through the equilateral trian- 
gle, $ = \ipi-2 + <^23 + <P3i\ = f • Thus we get non-zero 
spin chiral order parameter along these loops. See the 
results in Fig.2.(b). 

In summary, we have predicted an emergent chiral 
spin liquid state base on the topological Hubbard model 
on honeycomb lattice with spin rotation symmetry. In 
the end, we address the relevant experimental realization 
and the way to be conformed by numerical approaches. 
In condensed matter physics, there is no such material 
with a Hamiltonian of the topological Hubbard model. 
However, such system may be simulated in optical lat- 
tice of cold atoms. In Ref. [23, 24], it is proposed that the 
(spinlcss) Haldanc model on honeycomb optical lattice 
can be realized in the cold atoms. When two-component 
fermions with repulsive interaction are put into such op- 
tical lattice, one can get an effective topological Hubbard 
model. It is easy to change the potential barrier by vary- 
ing the laser intensities to tune the Hamiltonian parame- 
ters including the hopping strength (f-term) and the par- 
ticle interaction ([/-term). On the other hand, one may 
check our prediction by quantum Monte Carlo (QMC) 
simulations including the global phase diagram, the topo- 
logical degeneracy, the spin chiral order parameter [4]. 
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